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We study the electric flux tubes that undertake color confinement in A/" = 2 supersymmetric 
Yang-Mills theories softly broken down to A/" = 1 by perturbing with the first two Casimir opera- 
tors. The relevant Abelian Higgs model is not the standard one due to the presence of an off-diagonal 
coupling among different magnetic (7(1) factors. We perform a preliminary study of this model at 
a qualitative level. BPS vortices are explicitely obtained for particular values of the soft breaking 
parameters. Generically however, even in the ultrastrong scaling limit, vortices are not critical but 
C ' live in a "hybrid" type 11 phase. Also, ratios among string tensions are seen to follow no simple 

^ , pattern. We examine the situation at the half Higgsed vacua and find evidence for solutions with 

the behaviour of superconducting strings. In some cases they are solutions to BPS equations. 
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I. INTRODUCTION 

o , 

. Certainly, one of the most beautiful ideas in the context of quantum chromodynamics (QCD) is the confinement 
mechanism envisaged by 't Hooft [Q and Mandelstam [|| through the condensation of light monopoles. In essence 
. it states that the QCD vacuum should behave as a dual superconductor where magnetic order takes place, and 
"V* ' electric flux tubes form thus producing color confinement. In the context of QCD it stands for a kind of descriptive 
Oh, scheme, as long as it is not known how magnetically charged quanta can arise and condense in the effective low 
^ • energy theory. In this respect, the idea of Abelian projection proposed by 't Hooft has received increasing support 
from numerical simulations on the lattice in the last few years Even in the continuum, recent work using a novel 
parametrization of QCD jj] , points in the direction of the above scenario for color confinement |^ . From the analytical 
side, the understanding of non-perturbative phenomena in four dimensional quantum field theory has been put several 
steps forward since Seiberg and Witten constructed an exact solution for the low energy dynamics of SU{2) Af — 2 
d ' supersymmetric Yang-Mills theory In particular, it was possible for them to show that the mechanism of color 
confinement devised by 't Hooft and Mandelstam takes place when supersymmetry is broken down to Af — I. These 
results were soon extended to the case of SU{N) 0-^. Furthermore, when Af = 2 supersymmetry is softly broken 
down to A/" = 0, the same mechanism has been shown to persist ]To[ |. 

In spite of the fact that these results are well known, not much attention has been paid to the actual solutions in 
the strong coupling limit corresponding to electric fiux lines that would undertake quark confinement. In Ref. 
it was shown that this sort of vortices should have a spectrum of string tensions that distinguishes among different 
factors in the magnetic U{1)^~^ theory arising in the infrared. The same result was found in the framework of the 
M-theory fivebrane version of QCD, also named MQCD The string tension of the iV — 1 electric flux tubes 

Tfc , fc = 1, TV — 1 , is given -up to a dimensionful factor which is different for each theory but independent of 
k-, by a dimensionless function fN{k) = sui{Trk/N). This function is somehow universal as long as the soft breaking 
perturbation is carried by a single Casimir operator . Even in that case, the problem of finding such solutions 

in the particular model that emerges in this context has not yet been addressed in detail, |^ probably due to the 



^ Except for SU{2), both at the maximal singularity of the Coulomb branch for pure gauge jl3|, and on the Higgs branch in 
the theory with massive fundamental matter [|l4[. 
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naive expectation that the effective theory consists of A'^ — 1 copies of the standard U{1) Abehan Higgs model. It was 
already shown by Douglas and Shenker that the magnetic U{1) factors of the infrared quantum theory describing the 
neighborhood of the monopole singularities are coupled The existence of these off-diagonal couplings, r°^, was 
confirmed in two different frameworks. First, they appear in the expression of the Donaldson-Wittcn functional for 
gauge group SU{N) Moreover, these couplings were shown to satisfy a stringent constraint coming from the 

Whitham hierarchy formulation of the Seiberg-Witten solution in Ref. where, in addition, a general ansatz for 
T°f is given. 

In this paper, we extend the work to the case of SU (N), Af = 2 supersymmetric Yang-Mills theory softly broken 
to Af ^ 1. The analysis is performed in a "peculiar" scaling limit (named "ultrastrong" in We show that, even 

in that limit, generically there are no BPS electric flux tubes. A perturbative analysis leads to the conclusion that 
the phase of dual superconductivity is of type II i.e. there is a short range repulsive force between different vortices. 
This fact supports the expectation that indeed electric flux lines are safely confined into stable flux tubes, a feature 
of the confinement mechanism which is not a priori granted. ^ It is worth mentioning in this respect that numerical 
simulations in lattice QCD seem to point out that the type of Abclian Higgs model behind the picture of dual 
superconductivity is a critical one between type I and type II [|l8| . 

The plan of the paper is as follows. The setup of the problem is given in Section II where some aspects of the low- 
energy dynamics of Af — 2 supersymmetric gauge theories softly broken down to A/" = 1 are reviewed. We emphasize 
the existence of non-vanishing couplings between the different U{1) factors -even at the maximal singularity-, which 
play an essential role in our results. In Section III, we show that the string tension of vortex-like configurations obeys 
a Bogomol'nyi bound in the ultrastrong scaling limit. However, there are no BPS electric vortices in the system unless 
the complex phases of the soft breaking parameters corresponding to different Casimir operators are aligned. Even 
in this case, we show that the string tensions of the resulting BPS vortices are governed by a dimensionless function 
/7v(fc) which is different from the one obtained in ]^,pT|, the latter being recovered as a particular limit of our system 
corresponding to a single quadratic Af = 1 perturbation. In Section IV, we focus for convenience on the group SU (3) 
and analyze the critical vortex solutions in certain simplified cases. We speculate about the full spectrum of such 
configurations. A perturbative analysis of the dynamics expected for nearly critical vortices is performed in sections 
V and VI by means of energetic arguments. This analysis reveals the existence of repulsive forces among vortices 
corresponding to different magnetic U{1) factors. Thereafter we refer to this phase as an "hybrid" type II phase. 
In Section VII, the half Higgsed vacua are considered. The similarities and differences with the model proposed by 
Witten to describe cosmic superconducting strings [T9| ] are discussed. We find solutions to the Bogomol'nyi equations 
with the behavior of superconducting strings. Finally, Section VIII is devoted to our conclusions and further remarks. 



II. INFRARED DYNAMICS AT MAXIMAL SINGULARITIES 



The quantum moduli space of vacua A4\ of SU{N), Af — 2 supersymmetric gauge theory has a singular locus 
given by hypersurfaces of complex codimension one which may intersect with each other Along each of these 
hypersurfaces, an extra massless degree of freedom -whose quantum numbers can be read off from the monodromy 
matrix corresponding to a closed path encircling the singularity-, must be included into the effective action. At the 
intersections, many states become simultaneously massless. Of special interest are those singularities where — 1, 
i.e. the maximum allowed number of mutually local states, become massless. They are accordingly called maximal 
singularities. |^ 

The addition of a microscopic superpotential breaks supersymmetry and leads to an A/" = 1 theory 

^ 1 

= E I ^'^ Tr*'^ . (1) 

fc=2 



^ It could happen, for example, that the electric vortices result to be unstable, and their core grows and smears in such a way 
that they do not lead to confinement of electric charges jl7j. 
^ In SU{S), for example, the singular locus is given by the complex curves 

4u^ - 27{v ± 2A^f = , 

where u = l/2{Tr(fP) and v — l/3(Tr0^) are the gauge invariant order parameters constructed out of the scalar field belonging 
to the J\f — 2 vector supermultiplet, and A is the quantum dynamical scale. Higher intersections of these curves lead to the 
so-called Z2 and Z3 singularities, given respectively by the points {u^ = 27A^,v — 0} and {u — 0,v^— 4A®}. 
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Notice that those contributions in with > 3 are non-renormahzable. However, this does not necessarily mean 
that they do not affect the low-energy dynamics. They could be dangerously irrelevant operators [^,^. We will not 
discuss these subtleties here, and shall restrict ourselves to the case of up to cubic perturbations. 



1 o 1 , 

Wj^^i = - M Tr$2 + - Tr$3 

2 3 



(2) 



This breaking is soft, in fact renormalizable. The continuum vacuum degeneracy is lifted except for a given set of 
points which depend on the actual values of the parameters /i and ly. Q 

Let us focus on the low-energy effective field theory near a maximal point that we choose, for simplicity, to be that 
with real quadratic Casimir, u = NA"^. This is a dual Af = 2 supersymmetric gauge system with gauge group U{1)^~^ 
which includes both chiral multiplets ^'f — (xf;^^)j as well as hypermultiplets Hi = {Mi, Mi) that correspond to 
the monopoles that become light in that patch of the moduli space. One can choose a homology basis for the cycles on 
the auxiliary curve such that each monopolc has a unit of charge with respect to each dual gauge field. The quantities 
xf. Mi and Mi are chiral J\f = 1 superfields, while Vi^ are Af = 1 vector superfields (and W^^ their corresponding 
superfield strengths). For completeness, we give also the Af — content of these superfields 



where the notation for fermionic, bosonic and auxiliary components is the standard one. Setting af = {(pf), the 
coordinates at the point of maximal singularity that we are focusing on are af = 0. The dominant piece of the Af — 2 
low energy effective Lagrangian is given in terms of a holomorphic function called the effective prepotential 



— Im 

47r 



dH ^^^^-t 



d^eiMje 



M,, 



Mi\ + Re d^e W{x",M, M) 



(3) 



The monopole fields have been "integrated in" in order to soak up the singularity of the effective action when 

_ /iD\ 



where Mi becomes massless. The effective superpotential at low energies is 

W'(X^, M, M) = V2M,x?m + M i^ix"") + y vix"") , 



(4) 



the last two terms being the effective contribution of the supersymmetry breaking superpotential (^. In fact, U and 
V are the Abelian superfields arising respectively from the quadratic and cubic Casimir operators in the low-energy 
theory. The vacuum expectation value of their lowest components, U and Y , are the holomorphic coordinates in A^a, 
([/) = u and (F) = V. 

Written in component fields, the bosonic sector of the system is described by the Lagrangian 



'-eff, B 



D 



1. 



-hi,D,D, + h,F:F, 



Fnii Fnii 



where the auxiliary fields are solved as 



F^ : 

whereas field strengths and covariant derivatives are given by 



ij w-rrrij \ \4'mj\ ) 

-V2(t)f*4>* 



Frrn Fnii 



(5) 

(6) 
(7) 



* For instance, in the case of SU{3), the theory has generically five Af — 1 vacua, three of which are the maximal Z2 points. 
In the limit — > 0, the remaining vacua approach the Z3 points M. 
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(F^.). = 9^(^?),-a.(A^), , (8) 

Concerning Cj in (||), it stands for 

= ^iU,{^^) + z.y,-(0^) ^ iQle^ft- , (10) 

where Uj and are the derivatives of U and V with respect to 

C/,(</,^) = 4°)A + ^ii;PV^)Ai-f , uf = -2jsin^, (11) 

p>i 

y,(0^) = + ^.j''V)A^-^ , = -2jsin2^, , (12) 

p>i 

while u'^\(f)^) and v^^\(j)^) are homogeneous polynomials in of degree p, so that Cj are regular functions in the 
vicinity of the maximal singularity. Finally, bij is (^ times) the imaginary part of the period matrix r^*, 

^« = ^ = i + < + ° (x) ■ (") 

where Aj = Asin6'j and 9j = jir/N. When expanding around the vacuum expectation value af — {(l>f), T^{(j)-^) 
yields the effective coupling constant matrix. The logarithmic singularity when af — corresponds to the perturbative 
running of the dual coupling constant up to the maximal point, displaying the asymptotic freedom of the dual 
description. The coupling flows to zero due to the fact that the quantum fluctuations of massless monopoles have 
been integrated out. This is fine as long as one is interested only in searching for vacuum solutions. Then M and 
Aft in (|)-(|1) stand for the zero modes of the monopole field (see the discussion in [T0|). Here, however, in order 
not to run into double counting of degrees of freedom we should introduce, on physical grounds, an infrared cut off 
for the monopole loop integrals. In each U{1) factor the natural energy scale is set by the soft breaking parameters 

af - ICfV^' with 

^ = nuf^^A + vvf^K^ = -2iA(/xsin&, + M sin 2^0 , (14) 
and the perturbative couplings of each monopole to its corresponding dual vector field. 



9m 




(15) 



show logarithmic variations among different U{1) factors. ^ 



Even in the close vicinity of the singularity, different magnetic U{\) factors are coupled through t°^^ [QiSl- Exactly 
at the singularity, i.e. at af = 0, the generic expression proposed in [ p^ for these off-diagonal couplings is 

= sin fc^i sin fc^j cos k0p cos k0q , (16) 



k—1 P^Q—^ 



where Tpg' is given by 



4q = COS Op -2 COS f?fc)2 - (1 - Spq) log(2 cos Op- 2 cos 0q)\ (17) 

k^p 



^ In other words, we are dealing here with a macroscopic (classical) theory of the Ginzburg-Landau type, and we should 
consider the coupling constant of the Mi and M/ classical fields to V/^: wave-particle duality connects go with the running 
coupling constant of the quantum theory through the formula hgoi = gmiaf ~ 1^^-"' l^/'^), the strong coupling limit becoming 
the classical limit for the magnetically charged quanta [E2|. 
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with Op = {p — 1/2)tt/N and p,q = 1,...,N. In the case of SU(3), for example, r°|^ = i/7rlog2 R| 16|. These 
interactions are also present in the effective potential obtained from the terms in square brackets of (U), 

+ b^\^''){V2(j)mAm, + Qic^''))iV2(bm,4>m, + Ql^""))* • (18) 

Notice that, being positive definite, the potential is either positive or zero. Given the expectation values of the 
complex scalars, 

= af (0,„^) = rui {4>,n,) = rhi , (19) 

A/" = 1 supersymmetric vacua are in one to one correspondence with zeroes of Vef / : 

V2m,m, = -C,(a^) , (20) 
rriiaf — rhiaf = , (21) 
|m,| = \m\ . (22) 

i = l,2,...A'^ — 1. From ( ^l| ) we learn that monopole condensation can only occur at hypersurfaces where = for 
some i. At the maximal singularity, every af vanishes, and it is clear from (|20|)-(|2^) that — 1 monopoles have a 
chance to condense. While soft breaking is parametrized by /i and i/, monopole condensation is controlled by Ci. If 
for some j we have aj' = and adjust cj*^"* — 0, the corresponding U{1) remains unbroken {nij — rhj = 0), and the 
vacuum is said to be "partially Higgsed" . Summarizing, the Higgs vacuum Ti. at the maximal point is given by 



H = {m„TO, /|TO,;|2 = |m,|2 = |cfVV2, m, = -e'^^rm*} , (23) 

withCf^lcf^le^'^r given in equatio n (p^ . Since the absolute phases of are not fixed, it has the topology of 
a torus of genus g = N — 1. Equation (g3|) shows that the scalar components of the monopole superfields condense 
in the vacua placed at the maximal points. Although the presence of condensation suggests that confinement indeed 
takes place, some further analysis is required before this can be definitively established. An important question to be 
answered is whether the collimation of the electric (or dual magnetic) flux lines is energetically favored or not. This 
is a dynamical issue that goes beyond the simple vacuum analysis. 



III. BOGOMOL'NYI BOUND IN THE ULTRASTRONG SCALING LIMIT 

The resulting effective theory we have arrived at, in the bosonic sector, is an Abelian (TV — l)-IIiggs model with 
coupled U{1) factors and a quite non-standard Higgs potential. The search for stable vortex solutions in the complete 
system is a hard problem. On general grounds, one should not expect to have BPS string solutions in spite of the fact 
that Af — I supersymmetry is enough, generically, to have BPS vortices in four dimensions ||2^,|2^. At least, this is the 
case of TV = 1 QCD, where the strings are conserved modulo N so they cannot carry an additive conserved quantity 
such as a central charge |25|| . There is a limit, however, in which the system simplifies and admits BPS vortices 
prf . It happens whenever the condensation parameters ( p^ are independent of cf)^ , something that corresponds to 
linear perturbations in the superpotential (Q), i.e. Fayet-Iliopoulos terms. This kind of terms together with properly 
normalized quartic potentials are known to lead to Abelian Higgs models that admit BPS vortices |^,^,^. Taking 
into account (^l])-(|l^), one should consider A ^ oo and small values of the soft breaking parameters /i — s- and 

^' — > 0, such that /iA and vA^ remain finite. In this "ultrastrong" limit, C'i{4>D) C'^"'' are constants, and one can 
easily check that setting af — {(f)f) = is a consistent constraint. One may then study the existence of extended 
solutions in the remaining fields. 

The (bosonic part of the) effective Lagrangian adopts the following form: 



'1 



2 



(24) 
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where fe^^' stands for the actual value of 6^ at the maximal singularity and D'^^'' , F^"' are obtained from (|^) by 
replacing bij with b\^\ It is now feasible to give an expression a la Bogomol'nyi for the energy per unit length 
corresponding to static and magnetically neutral = 0) vortex-like configurations (i.e. configurations with 

translational symmetry along one axis) by means of the remainder Af = \ supersymmetry [£3|,^ (see also where 
a multi Higgs system has been treated). Indeed, the energy density can be rearranged as follows: 



(0) p(0) 



2 'J V 



FI2 ± 



(0) 



(Fi 



(0) 



(25) 



where the last term, corresponding to the current Jb = —i{(j)'^.Db(f>mi + i'm ^bi'mi), does not contribute to the string 
tension for finite energy configurations. It is easier to analyze this system in a different set of variables, obtained from 
the above ones by means of an SU{2)fi transformation yielding 



V2F^°^ 



Df^ =-V2Re(e'"^;^°^) , 



V2' 



+ e" 



(26) 

(27) 
(28) 

(29) 



The tension (Jeff = J ^eff now reads 



aeff = / d^x [\bf [fi, ± A^")) {fi, ± i)f ) + \{D, ± ^D: 



2)(Pmi 



(0)n 



(30) 



The last term breaks explicitely SU{2)ji symmetry. Finiteness of the string tension demands regularity of the fields 
on M^, and vanishing of the potential energy, field strenghts and covariant derivatives at infinity. Altogether, these 
requirements make the space of solutions to split into Z^^^ disconnected pieces that differ by the winding numbers of 
each (prui over the border of the plane. The electric-fiuxes label these sectors. In particular, in the (ni, 712, titv-i)- 
sector they are 



F(, = 2iTn 



1 



j = l,2,...,7V_l 



(31) 



The string tension of possible vortex configurations with topologically quantized (ni,n2) electric fiux, exhibits a 
Bogomol'nyi bound 



o-e/ / > 4V27rA ^ I (^ sin 9, + vK sin 29,) cos(a + (if' ) 



(32) 



which is saturated for configurations solving the following set of first order equations: 



FI2 = ±V2 Re(e^"Ff ^) , 
[Di ± iD2)im, = , 
The second equation in ( ^ implies 

|0mj = |0mj| Im 



A"' +iV2 Im(e*"i^^°') = , 

(i?l±ii?2)L. = . 



(0) 



= . 



These constraints should hold at any point, in particular, at zeroes of the Higgs field. Thus 



(33) 
(34) 

(35) 



(36) 
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with a + j3- — or tt. Consequently, for i j, = (3l — (3 - —Q or tt. Summarizing, there are no BPS electric 

vortices in the system unless the complex numbers cf'^ are aligned or anti-aligned. This aUgnment, in turn, requires 
supersymmetry breaking parameters to have no relative complex phases. Notice that this corresponds to having a 
CP invariant bare Lagrangian. For definiteness, in the case of S'[/(3) one easily sees that 

Cf ^ = %/3A(^ + uK) ■ ^ = V3A(/i - vK) , (37) 

so that /Sji'' = or TT if and only if aig[vA) = arg/i + n-K and \vh\ < \fi\ or \i/A\ > \fj,\ respectively. 

A comment is in order at this point regarding the string tensions of unit vortices, whose existence will be discussed 
below. It is immediate to read, from (^), the string tension of electric vortices carrying a single flux quantum 
rifc = 1, n^^fc = 0. Up to a common factor, it is given by 

Tk cx A/w(fc) /jv(fc) = \fi sin 4 + i^A sin2^fc| . (38) 

This result makes clear the dependence of ■/ 'jv( fc) on the supersymmetry breaking deformation entering the superpo- 
tential. It generalizes previous results in ||q, |ll| , [l2| and, in particular, it shows that for perturbations other than the 
quadratic one, the string tensions arc modified with respect to those in the above mentioned results. In particular, 
notice that when /i and v do not vanish it is possible to have different string tensions even in the case of SU{3) and, 
in general, Tj, 7^ Tjv.fc. 



IV. ALIGNED VACUA: CRITICAL VORTICES 



We will focus hereafter on the case of SU{3). When the constants /i and ly are fine tuned in such a way that the 

phases of the two complex energy scales Cj;*^-* and C2^^ are either aligned or antialigned, i.e. (3^^ = or tt respectively, 
we are at the self dual point. The Bogomol'nyi equations (^-(|3^), after (^, read 

Fi2 = ±lb[f-'e,i\ip,\'-v]) , (39) 
(Di ± itjD2)ifj = , (40) 

where = e^^^+'^f ^ = ±1 and bf^ is 

Also, we have performed, for convenience, some redefinitions of the fields, ipj — 2(f)mj , and parameters, w| — 2\/2|Cj"'' |. 
Let us further remark that Eq.(|39|) gives an unusual contribution to the electric field of each dual U{1) factor from 
zeroes of both Higgs fields. This is a straight consequence of the presence of off-diagonal couplings and leads to 



interesting results. It is clear that solutions to (39)-(|40|) also satisfy the Euler-Lagrange equations. Without loss of 

■MO) 

generality, we can adjust a so that ei = +1, €2 = e = e^'^^i = ±1. Let us focus on the BPS solutions with upper sign. 
The first order system can be written as 

^1^2 = \ii\ipi\^~vf)~e^{\ip2\^-vl) , (42) 

^f2--7(l^i|'-^^?)+eA2(|^2p-«|) , (43) 

(Di + iD2) - , (44) 

[Di + ieD2)^2 = , (45) 

with 

7 = &ir' = ^(3l.,iA2 + <?K2Ai). (47) 
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Note that, as we are in the weak (715-coupHng regime, 7 < Xi. Naively one would suspect that in the scaling limit we 
are interested, the system diagonalizes. Notice however the important fact that the relative factor between Xi and 7 
vanishes only logarithmically. Hence, for example, setting |Ci|/A? ~ 10^^° in ( p^ yields 7 ~ (log2/5)Ai ~ 0.13 A^. 

The topology of the configuration space determines global properties of the solutions in two ways: the quantization 
of the fluxes is due cither to the asymptotics of the Aj fields or to the existence of a prescribed number of zeroes of the 
ifj. These global inputs should be made compatible with the differential equations, as it happens in the Abclian Higgs 
model. In the present situation things are less clear; from Eqs.(^^-(^5|), where no mixing between both U (l)s shows 
up, one reads the electric fluxes using Stokes theorem and the asymptotics of Aj. On the other hand, Eqs.(^2|)-(|4^) 
mix the factors and both (fi and (p2 contribute together to each F^2. In this respect, our system is quite awkward as 
compared with other non-diagonal models as, for example, non-relativistic non-abelian Chern-Simons theories [p9| , 
in which the same mixing appears in the field strength and covariant derivative equations. Here, there is mixing in 
the former but not in the latter, and given such an asymmetry, it is much more difficult to show whether the local 
equations and the global conditions reconcile or not. 

On general grounds, it is reasonable to expect that the equations p^)--(^) will exhibit solutions in the topological 
sector (ni,rt2) with ni,n2 representing the integrated flux of an "ensemble" of noninteracting vortices located at 
different (maybe coincident) positions. Indeed, the smallness of the ratio 7/Ai suggests to consider this system as a 
pertubation of the diagonal situation, so that the above solutions would come out from continuous deformations of 
the standard critical Abrikosov vortices. Only in some simple cases, the question about the existence of solutions can 
be answered by taking advantage of known results from the standard Abelian Higgs model. This will be done in the 
following two situations 

• Solutions of type (rt, 0) and (0,n). 

Clearly it will be enough to prove existence of one type, say (n,0). Assume therefore that (p2 = |<P2|e'^^ is 
nowhere vanishing on the finite transverse plane. As usual, (E5[) couples ^2 and A2. So, if \(p2\ has nowhere 
a zero, regularity of the phase enforces A2 to have vanishing circulation around any loop. By Stokes theorem 
= everywhere, an inserting this back into ( p3| ) yields a constraint that correlates the profiles of \ipi\ and 

\V>2\, 

l¥'2p=ef + (48) 

A2 

Existence of the required vortex profile for can be proved by inserting ( ^ ) into (|4^), which leads to the 
standard Bogomol'nyi equations for the critical Abelian Higgs model (after a suitable re-normalization of the 
Higgs field) 

^^i'2-Ai(^l-Q(|^ip-«?), (49) 
{Di + iD2)ipi = . (50) 

We learn from ( ^ ) that if \(pi\'^ ranges from (at the origin) up to (at infinity), |(^2p will correspondingly 
interpolate between —e-^vf + u| and To remain consistent with our initial asumption that \ip2\ vanished 

nowhere we must set either (3^^ — with w| > -^vf, or else e = —1, i.e. , = tt. We observe that the latter 
possibility is less contrived. 

• Solutions of type (n, n) for a single perturbation. 

Let us briefly consider the case of SU{3) Af — 2 supersymmetric Yang-Mills theory softly broken to A/" = 1 only 
by means of a single Casimir operator, i.e. fi — or :^ ~ 0. In both cases, = or tt, and the theory is 

critical. Moreover, Ai = A2 = A, cj"-* = c!f ^ and hence vi = V2 = v, so that the Bogomol'nyi equations have an 
almost trivial solution of vorticity (n, n) (or (n, —n)), by imposing the ansatz (pj = Lp, Aj = A (or = V'l = V'j 

— A2 = Ai = A) in the case (3^^^ = (or tt) . The system is again reduced, after a suitable normalization of the 
Higgs field, to the critical Abelian Higgs model 

^^i2-±(A-e7)(|¥>p-z^') , (51) 
[Di ± iD2)Lp = , e = e'^'^i' . (52) 

It is crucial, for the system to admit regular solutions, that 7 < A as it indeed happens. As it is well known, 
the general solution to this sytem represents an assembly of n separated vortices centered at the zeroes of ip. In 
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our case, every such zero is doubled and we have assembhes of n couples of superimposed vortices of both U{\) 
fields. 

Also, self-dual configurations in which the center of the vortices of different types split apart, can be easily 
constructed along the lines in |30|j3l| ]. To see this, we perturb one of the solutions just described for (3^^ — 

^p'.j = tpj + Sipj , A'j = Aj + SAj , (53) 

and linearize the self-duality equations to get 

- AidJAi - 2Xip*Sipi + 2jip*6ip2 = , (54) 
-4idJA2 + 2jip*Sipi - 2Xip*6ip2 = , (55) 
igofSAj + {dg + igDAj)Sipj = , (56) 

where we use the notation dz — ■^{di — 182), Aj = ^[{Ai)j + i{A2)j], j = 1, 2, and fix the gauge conditions as 

dc{6A,)i = -Xlifl^Sn^ + j\ip\''dn2 , (57) 
dc{SA,)2 = M^l'^Sfli - \\ipfdn2 . (58) 

By writing 6ipj — (p^j and using (|56|), the vector perturbations are found to be SAj = j^^z^j and the system of 
linearized equations reduces to 

V^W± ^ 2{Xt i)9d\v^\^W± , (59) 

with W± = ^1 ± ^2- Notice that in both equations (A ^ 7)50 > 0. Although they have not regular square- 
integrable solutions, we can admit singular ones provided the singularities of fit with the zeroes of (p in such a 
way that Sipj is well-behaved. Take for instance the case of a radially symmetric solution of vorticity n centered 
at the origin of the complex plane. Then, for small z 

ip{z,z)^z^, (60) 



and a singularity of W± at the origin is harmless if its order is lower or equal than n. Equation (|59|) has indeed 
solutions with such a behaviour | |3^ . To be exact, two sets of linearly independent self-dual perturbations 

W^{z, z), m = 1, 2, 3, , n with 

W^iz,z)~z-"\ z~0. (61) 
In particular, if we consider W± = —aW^, we get, near the origin, 

6 ^ -az-"" 6 - , (62) 



so that 



ip[ ~ z"-"(.2™ - a) ip'2 ~ .2" . (63) 



This perturbation realizes the splitting of a (n, n) vortex at the origin into a (n — m, n) at that point and m 

(1,0) vortices located at the m roots of the coefficient a. The analysis for P^^^ = n (i.e. e = —1) is totally 
equivalent and yields nothing but vortices of type 1 and anti- vortices of type 2 or viceversa, moving freely with 
respect to each other. 

For the general analysis, following Jaffe and Taubes |Q, the Higgs fields should be "couched" as 

= Vje5(«^+'"^) , (64) 

to recast the Higgs system in the following form 

V^ui = 2Aiz;2(e"i - 1) - 2e7'y2(e"^ - 1) + Ebc^c'cf^i , (65) 
V^U2 = -27'y2(e"i - 1) + 2eA2'y^(e"^ - 1) + ebcdbdcn2 . (66) 
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The gauge fields are determined by 

{A,)i^-'^{d,n^+e,adaui) , (67) 

{Ac)2 = -^^{dc^2 + ecadaU2) ■ (68) 

At each (ni, 77.2) sector, regularity impHes that ipj has exactly rij zeroes on C, say z^J, . . . , z^ . . Also, these are the 
only points at which the singularities of the phases can occur. We can then choose the particular gauge 

rij 

^■(z,z) = 2^arg(z-z/) , (69) 

in which the problem reduces to 

V=^ui = 2Aiw?(e"i - 1) -2e7i;|(e"^ - l) + 47r^(5(z-z/) , (70) 

1=1 

V^ua = -27t-2(e"i - 1) + 2eA2W2 (e"^ - 1) + 47r ^ 6{z - zf) , (71) 

where both Uj should vanish at space infinity. The general analysis is involved, and usually goes through by numerical 
relaxation techniques or hard Sovolev estimates. 

V. HYBRID TYPE II VORTICES 

By itself, the Abelian Higgs model we are dealing with is worth a detailed analysis. For the moment, and awaiting 
a sounder analytical or numerical study of its solutions, aside from the two simplified samples considered above little 
can be said about the generic (ni,n2) vortex solution. An interesting peculiarity comes from the fact that there are 



only two overall choices of signs available in equations (p9) and (40): either upper or lower sign have to be taken 
simultaneously on all the equations or, else, the bound (^) will not be saturated. This should be contrasted with 
the situation in the standard diagonal Abelian Higgs mode , where each C/(l) can be conjugated independently. To 
better grasp what is going on let us consider the Bogomol'nyi equations (^2[)-(p5|) with f]^^^ = 

^ ±{X,W, - , (72) 

{D^±iD2)^i^0 , (73) 

Fi^^±{X2W2-jWi) , (74) 
(Di ± iD2)^2 = , (75) 

with Wi — (IViP ~ ■ 7 '^ii-^2j (i"-iji'^2) vortex with ni,n2 > come from solutions to the previous 
equations with the upper (lower) sign which should correspond to deformations of analogous configurations in the 
case 7 = 0. In the diagonal limit 7 = the vortcx-antivortex solutions (±ni,=pr7,2) would also solve the previous 
equations but with a choice of sign for (|7^)-([73|) and the opposite one for (]7^)-([75|). If 7 ^ 0, as is now the case, 
solutions with this second choice of sign do not saturate the bound (^) and, indeed, there is an energy remnant 
coming from the off-diagonal piece £ — tt \nivf + riai'il + 

For anti-aligned magnetic fields, this extra term is negative and tends to increase the overlap by attracting the cores 
of vortices of different kind. 

A similar reasoning can be carried out of p'^i = tt. In this case, the equations read 

=±(Ai 1^1+ 71^2) , (77) 
{Di±iD2)vi^Q , (78) 

F^2^T{\2W2+lWi) , (79) 
{Di T iD2)^2 = , (80) 
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and critical configurations are naturally of the form (±ni, ^"2), n-i,n2 > saturating the bound £ = tt \nivf — n2yo \- 
Here, in contrast, vortex-vortex solutions of the form (±ni,±n2) would lead to the same energy surplus as in (|7q). 
But now 5£ > for aligned magnetic fields, and this term decreases by minimizing the overlap, hence by taking the 
cores far appart. 

In summary, to a first approximation, we see that, if not neutral, vortex-vortex (vortex-antivortex) configurations 
behave repulsively (attractively) as in type II superconductors. Since this interaction involves vortices of different 
C/(l)'s, we speak of an "hybrid type 11" phase. 

Let us discuss the peculiarities that arise whenever one tries to model confinement in the present scenario. First we 
fix some notation for convenience: the chromoelectric fluxes (^1,^2) of the basic vortices arising in the dual Meissner 
effect are (1,0) ("vortex 1") and (0,1) ("vortex 2"). In turn, quarks enter the system as external probes with 
chromoelectric charges {Qi,Q2) equal to (1,0) ("red quark"), (0,-1) ("blue quark") and (—1,1) ("yellow quark"). 
{hi, /i2) is the "monopole" basis of the Cartan algebra of the dual SU{3) group and the fundamental BPS monopoles 
correspond to the simple co-roots of SU{3). In other words, the chromomagnetic charges of the fi-iield quanta is 
hi ~ l,hj^i — 0. Consider now, for example, the case f3^i = 0. According to our previous analysis, chromoelectric 
fiux tubes of both (1, 0) and (0, 1) type form in response to parallel external electric fields Ei and i?2. Vortices of type 
1 end at pairs of red quark-antiquark and vortices of type 2 finish at pairs of blue antiquark-quark. There is therefore 
confinement of red and blue quarks in a critical phase between Type I and Type II superconductivity, whereas the 
yellow quark confinement occurs in a hybrid Type II phase. The weak repulsion between the vortex 1/antivortex 2 
pair pull slightly apart the fiux lines from each other. Thus, the quark/antiquark potential energy would increase 
slower than linearly with the distance, and one is allowed to expect deviations from the area law, but the force is still 
confining. If, instead, (32^ — tt, a pair of yellow quark-antiquark will now be joined by a stable and non-interacting 
vortex 1/antivortex 2 pair of flux tubes. In conclusion, the cases = or tt can be physically distinguished by 
the behaviour of the yellow quark-antiquark force. At large separation W-pair production leads to instability of the 
string and the lowest string tension governs the large distance regime [p|JTl|] . 

In the framework of condensed matter it is well known the fact that, in standard type II superconductivity on a 
finite piece of material, though mutually repelling, vortices tend to form a regular pattern by lying at the sites of a 
triangular lattice. This fact can be reproduced analytically by variational methods p^ . We expect a similar situation 
here, the difference being that now repulsion involves vortex cores of distinct Higgs fields. Upon substitution of ( |3^ ) 

into (|2^), the exact second order equations with f]^^^ — tt, corresponding to vortices of type 1 and 2, in a finite piece 
of material 



bf[^da{Fab)l+b^°2^da{Fabh = 




- LpiDblfil) , 


(81) 


bf^da{Fab)2+b^2lda{Fab)l = 


^{Lp^Db(p2 


- <^2-DfcV2) : 


(82) 




1 *J,(0) 




(83) 


DcDcip2 = 


1 *l(0) 
-^^^2^2, 




(84) 



should now be supplemented with periodic boundary conditions. Thus, the system of differential equations is defined 
in a torus of modular parameter t = L2/ Li e*^. We have chosen the xi-axis as the direction of the first Li periodicity; 
the length and direction of the second periodicity is determined by ^26*^. Application of the Rayleigh-Ritz variational 
method as in plus previous work on the role of Riemann Theta functions in magnetic systems [ p5[ , suggest the 
field configurations 



ipi = C„ij exp[mimilmz — -(Rez — nimi)^] , (85) 

miGZ 

C„i2 exp[m2»Ti2lmz — -(Rez — 712^2)^] , (86) 



V52 = 

m2eZ 

where ni , 77.2 are integers and z = y^g]j{X — 7) ( — - j , as trial functions to model extremals of the energy. In 

fact, the choice of the coefficients Cmi and in such a way that 
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ip^'^{z) — exp{— 7rni 
fTi^) = exp{-7m2 



(Im^)2 ^ 
-} 

h=i 



Imr 

(Im^)^ 
Imr 





ii. 



(^1-) , 

ni 



(2=1 



, 2 

is. I 1 

"2^2 



(^1-) 

"2 



(87) 
(88) 



leads to (meta)-stable solutions to the field equations. Here — l,...,ni, and Oi^K^^lr) are the Riemann Theta 
functions with characteristics, see [p5[ and references quoted therein. 
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FIG. 1. The Type II hybrid lattice. Black and white circles represent the core of vortices corresponding to different U{1) 

Notice that the solution describes ni chromoelectric vortices, located at the zeroes of ip^^ , and n2 vortices of the 
other kind centered around the zeroes of (p2^ . It corresponds therefore to a hybrid static triangular lattice of vortices; 
see the Figure. One can check, from a dynamical point of view that a configuration like this where a vortex of type 1 
is at the center of a square with vortices of type 2 at the vertices and viceversa is stable against small fluctuations. 



VI. MISALIGNED VACUA 



As discussed earlier, there are no BPS vortices in the generic case where soft breaking parameters are not aligned. 
We would however be interested in the response of the BPS configuration when an infinitesimal misalignment f3^^ = e 
or = TT + e is turned on. The dynamics of the system drives the configuration off the constraint ( |36| ) which, 
therefore, can no longer be imposed consistently. In fact, though the Higgs mechanism yields a critical mass spectra 
for any value of e (an obvious consequence of supersymmetry), the eigenvectors do depend on this phase difference in 
such a way that when it is different from or tt, massive excitations do not respect the constraint surface (^6|). 

In the same vein as above, for small values of e we will treat the system as a perturbation of the critical situation 
in which the net effect of the misalignment reflects itself in a force between the former noninteracting vortices. The 
shortcut to obtain the sign of this force is to split the energy ( pOf) of the configuration as a BPS contribution plus 
an additional perturbation. Namely, after inserting the ansatz ( p6| ) into (24), solutions to (^9|)-(40) exhibit a string 
tension aef / — cTeff + ^'^ef / where cr^jj is given in ( ^2|) and 

6a,ff^e^Jd^x{\ip,\'-vl)i\^2\^-vl). (e « I) (89) 

with e = 1 for f3^ = + e and e = —1 for /32i'' = ir + e. Consider a vortex configuration of type (1, 1) where the zeroes 
of each Higgs field are well separated. Then, the above surplus of energy is positive for e = 1 and decreases as the cores 
are taken further appart and the overlap diminishes hence the interaction in this case is repulsive. When perturbing 
around the anti-aligned case, P^2i = 7r + e, the energy increment ( ^9| ) reverses sign. Previously non-interacting, (1, —I) 
antiparallel vortex configurations tend to increase the overlap in order to lower the perturbation, and hence the force 
is attractive. 
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In summary, when perturbing around the ahgned or misahgned scenarios, the vortex configurations are not neutral 
anymore, and the interactions follow the pattern that was previously named "hybrid type 11" where, if made of distinct 
U(l)'s, parallel vortices repel and antiparallel vortices attract. 



VII. HALF HIGGSED VACUA 



As pointed out in j^, for particular values of the soft breaking parameters and v we have four instead of five 
vacua. This happens whenever one of the two half Higgsed vacua {af* Ojof' = 0} with Ci(/i, j^) = 0, or (1 2), 
meets and replaces one of the normal vacua at {of = 0, = 0}. This possibility is actually achieved by turning off 

cf^ for i = 1 or 2. Since precisely at the Z2 point we have (|37|), this amounts to /i and v fulfilling /i = =fj/A. Let us 



choose for definiteness, cf'' = 0. Inserting this back into (18), the effective potential at the maximal point reads 



V = iAid^il^ - vlf + ^XM' - i7Cos/3i|^2|'(|^i|' - vl) . (90) 

Observe that the phase of ip2 is free. When cos (3i < this is precisely the type of situation that was studied by 
Witten ||l^ and shown to lead to superconducting strings for specific ranges of parameters. Let us briefly recall the 
essence of the mechanism. As the vacuum equations ( pO| ) exhibit, only the first U{1) is broken by the v.e.v. (ipi) — vi, 
whereas the second U{1) remains intact since ((/?2) = 0. This is fine for vacuum solutions, but suppose now that ipi 
developes a vortex line. At the core of the vortex {(pi) — and, in turn, it may become favorable that {ip2) ^ there. 
Actually the model considered in |l^ is slightly more general than ours involving the potential 

V ^ iffd^iP - v'f + ^.9l^2|^ + fW\^2\' - m^\V2\' ■ (91) 

The detailed analysis of the dynamics showed that for parameters in the range /w^ > m^, instability actually takes 
over and the superconducting string indeed forms. We see easily that the present situation lies precisely at the 
boundary of the region of validity, since in our case fv^ — rri^ = 0, and the induced mass term for (p2 exactly vanishes. 
In , this situation was also studied and seen to yield a power law decay of the profile of (p2 which leads to a long 
range scalar attractive interaction among vortices. 

At this point we would not like to put forward too strong a claim, but simply point out the ocurrence of this 
coincidence among models. The possible existence and relevance of structures like superconducting strings in the 
microscopic context of confinement models should be handled with care. For example the question of quantum 
tunnelling will be certainly much more relevant here than for cosmic strings. Incidentally this question was also 
addressed in p6| where it was seen that these power law solutions are more stable than the usual ones. 

As compared with Witten's model, the one here involves the additional feature of the non-diagonal kinetic term 



for the (dual) vector particles (c/. eq. (24)). But precisely the fact that the quadratic forms of kinetic term and 
potential are related paves the way to the possibility of rewritting the energy as a sum of squares (|30|). We may 
therefore expect vortex solutions of the superconducting type with dynamical properties of BPS configurations. We 
can check that this is indeed the case by looking at the smooth deformation of a generic (anti-)aligned scenario. ^ 

Let us follow a continuous line of anti-aligned {P^i — tt) vacua cf'^ ^ 0, clf'^ — + 0. Precisely in this situation, ( |4^ ) 
presents no obstruction to a smooth deformation of the (n, 0) solutions down to the situation where V2 — 0. In this 
limit the profiles of \(pi\ and \(p2\ are correlated in such a way that both vanish at opposite ends. In fact, as {ipil"^ 
varies from zero up to vf far away, \(p2\ interpolates between j^vf = (log 2/87r^)g|) ^vf at the origin (which need not 
be small !), and at infinity. Moreover, since the phase of ip2 is free, the same arguments of ref. |l^ can be used to 
show that a persistent current occurs. We would call this a BPS superconducting string solution. 



VIII. CONCLUDING REMARKS 



The present paper is devoted to the low energy dynamics of Af — 2 supersymmetric gauge theories softly broken 
to = 1 by a superpotential containing up to cubic perturbations. The effective lagrangian in the neighborhood 



^ As we approach the situation when C2 — » 0, the parameters that enter (^0|) are such that 7, A2 << Ai (see eqns.(|l5|) and 
(^6[)-(|47|)). Hence at very low energy the second U{1) seemingly decouples. This is suggested by the N — 2 exact effective 
solution, although it is reasonable to expect modifications of the renormalization group fiow in the TV = 1 theory. 
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of maximal singularities of the quantum moduli space corresponds to an Abelian U{1)^~^ multi Higgs system with 
couplings among different dual U (1) factors. The case of SU{3) has been analized in some detail. There are generically 
no BPS electric vortices in the system unless the soft breaking parameters have coincident complex phases (or they 
differ by tt) and the ultrastrong scaling limit is taken. We have seen that the effect over a BPS configuration of 
turning on an infinitesimal misalignment among these parameters is the appearance of a net repulsive force between 
parallel vortices corresponding to (zeroes of) different Higgs fields. In a finite piece of material, metastable solutions 
take place and vortices develope static triangular lattice. We call this phase "hybrid Type U" dual superconductivity. 

When the theory is perturbed with a cubic superpotential, the ratio of string tensions differs from that computed 
in the quadratic case Q both when the Tr<I>^ perturbation is present or not. In the former case, we found that 
these ratios even depend on the supersymmetry breaking parameters. These results were obtained after imposing the 
ultrastrong scaling limit. It would be certainly interesting to know if similar results emerge in the context of MQCD. 
This is intriguing in the sense that string tensions in MQCD are given by the distance of D4-branes which, for a 
single Casimir perturbation, are stretched at the roots of unity over a circle of radius of order A pT| , so one would 
not expect them to be modified (except, possibly, for a global factor due to an induced change in vVjas compared to 
the purely quadratic case. 

A natural extension of the present work involves the case of A/" = 2 supersynimetric theories softly broken down to 
J\f — 0, and possible soft breaking by higher than the two first Casimir operators. This program can be addressed 
within the Whitham approach to the Seiberg-Witten solution, where the slow-times of the hierarchy can be used as 
spurionic sources of soft supersymmetry breaking ||2l[| . 



ACKNOWLEDGMENTS 



We are pleased to thank Jose F. Barbon, A. Gonzalez- Arroyo, Michael Douglas, Amihay Hanany and Marcos Marino 
for interesting discussions. J.M. wants to thank J.J. Blanco Pillado for pointing out reference |^^. The work of J.D.E. 
has been supported by the National Research Council (CONICET) of Argentina and the Ministry of Education and 
Culture of Spain. The work of J.M. was partially supported by DGCIYT under contract PB96-0960 and European 
Union TMR grant ERBFM-RXCT960012. 



[1] G. 't Hooft, in High Energy Physics, Proceedings of the European Physical Society International Conference, Palermo, Italy, 

23-28 July 1975, Ed. A. Zichichi, Editrice Compositori, Bologna, 1976. 
[2] S. Mandelstam, Phys. Rep. 23 C, 245 (1976). 

[3] See, for example, M.I. Polikarpov, Nucl. Phys. Proc. Suppl. 53, 134 (1997) and references therein. 
[4] L. Faddeev and A.J. Niemi, Phys. Rev. Lett. 82, 1624 (1999); Phys. Lett. B 449, 214 (1999). 
[5] Y.M. Cho and D.G. Pak, Magnetic confinement in QCD, hep-th/9906198. 

[6] N. Seiberg and E. Witten, Nucl. Phys. B 426, 19 (1994) [Erratum ibid. B 430 (1994) 485]; Nucl. Phys. B 431, 484 (1994). 

[7] A. Klemm, W. Lerche, S. Theisen and S. Yankielowicz, Phys. Lett. B 344, 169 (1995); P.C. Argyres and A.E. Faraggi, 
Phys. Rev. Lett. 74, 3931 (1995). 

[8] M.R. Douglas and S.H. Shenker, Nucl. Phys. B 447, 271 (1995). 

[9] P.C. Argyres and M.R. Douglas, Nucl. Phys. B 448, 93 (1995). 
[10] L. Alvarez-Gaume, J. Distler, C. Kounnas and M. Marino, Int. J. Mod. Phys. A 11, 4745 (1996); L.Alvarez-Gaume and 

M.Marino, Int. J. Mod. Phys. A 12, 975 (1997). 
[11] A. Hanany, M.J. Strassler and A. Zaffaroni, Nucl. Phys. B 513, 87 (1998). 
[12] M.J. Strassler, Prog. Theor. Phys. Suppl. 131, 439 (1998). 
[13] W. Garcia Fuertes and J. Mateos Guilarte, Phys. Lett. B 437, 82 (1998). 
[14] A. Yung, Nucl. Phys. B 562, 191 (1999). 

[15] M. Marino and G. Moore, Comm. Math. Phys. 199, 25 (1998). 

[16] J.D. Edelstein and J. Mas, Phys. Lett. B 452, 69 (1999). 

[17] A. Achucarro, M. de Roo and L. Huiszoon, Phys. Lett. B 424, 288 (1998). 

[18] See, for example, G.S. Bali, C. Schlichter and K. Schilling, Prog. Theor. Phys. Suppl. 131, 645 (1998); Nucl. Phys. Proc. 

Suppl. 73, 638 (1999) and references therein. 
[19] E. Witten, Nucl. Phys. B 249, 557 (1985). 

[20] D. Kutasov, A. Schwimmer and N. Seiberg, Nucl. Phys. B 459, 455 (1996). 



14 



[21] J.D. Edelstein, M. Marifio and J. Mas, Nud. Phys. B 541, 671 (1999). 

[22] See p. 206. in S. Coleman, Aspects of Symmetry, Cambridge University Press, 1988. 

[23] Z. Hlousek and D. Spector, Nud. Phys. B 370, 143 (1992); Phys. Lett. B 283, 75 (1992); Nud. Phys. B 397, 173 (1993). 

[24] J.D. Edelstein, C. Nunez and F.A. Schaposnik, Phys. Lett. B 329, 39 (1994). 

[25] E. Witten, Nud. Phys. B 507, 658 (1997). 

[26] J.D. Edelstein and C. Niificz, Phys. Rev. D 55, 3811 (1997). 

[27] J.D. Edelstein, C. Nuficz and F.A. Schaposnik, Nud. Phys. B 458, 165 (1996). 

[28] E.B. Bogomol'nyi, Sov. Jour. Nud. Phys. 24, 449 (1976). 

[29] G. Dunne, Self-Dual Chem-Simons Theories, Lecture Notes in Physics Vol.36, Springer 1995. 

[30] E. J. Weinberg, Phys. Rev. D 19, 3008 (1979). 

[31] P. J. Ruback, Nud. Phys. B 296, 669 (1988). 

[32] T.M. Samols, Comm. Math. Phys. 145, 149 (1992). 

[33] A. Jaffe and C. Taubes, Vortices and Monopoles, Progress in Physics, Birkhauser, 1980. 
[34] A. Abrikosov, Fundamentals of the theory of metals, North Holland, Amsterdam, 1988. 
[35] J.M. Guilarte, J.M. Porras and M. de la Torre, J. Geom. Phys. 27, 297 (1998). 
[36] D. Haws, M. Hindmarsh and N. Turok, Phys. Lett. B 209, 255 (1988). 



15 



